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ABSTRACT 

We have created synthetic neutral hydrogen (Hi) Galactic Plane Survey data cubes 
covering 90° ^ I ^ 180°, using a model spiral galaxy from SPH simulations and the 
radiative transfer code TORUS. The density, temperature and other physical parame- 
ters are fed from the SPH simulation into TORUS, where the Hi emissivity and opacity 
are calculated before the 21-cm line emission profile is determined. Our main focus 
is the observation of Outer Galaxy 'Perseus Arm' H I, with a view to tracing atomic 
gas as it encounters shock motions as it enters a spiral arm interface, an early step 
in the formation of molecular clouds. The observation of H I self-absorption features 
at these shock sites (in both real observations and our synthetic data) allows us to 
investigate further the connection between cold atomic gas and the onset of molecular 
cloud formation. 
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1 INTRODUCTION 

The formation of molecular clouds out of cold, dense atomic 
gas is a key process in the cyclic evolution of the interstel- 
lar medium. Above column densities of Nh ~ W 20 cm~ 2 , 
atomic hydrogen gas transitions to molecular H2, a critical 
step towards the formation of stellar systems (|Dixon et al.l 
1 19981 ). For the molecular fraction of the gas to remain stable 
depends largely on the ability of the molecular gas to self- 
shield against dissociating photons from the local radiation 
field. Moreover, the importance of dust grains in H2 forma- 
tion requires their presence in these transition interfaces. 

The 21-cm neutral hydrogen line is the main tracer 
of interstellar structure. The ubiquitous nature of atomic 
hydro gen in the Galax y was demonstrated in early sur- 
veys (|Oort et all Il958l ). and more recent surveys imag- 
ing the HI distribution at moderately high resolu- 
tion (parsec scales at distances of kpc) have steadily 
improved upon that view. Almost all of the HI in 
the Galaxy's disc has been sampled by the Interna- 



vastly on the previous generation's all-sky HI maps. Galac- 
tic 21 cm surveys at Arecibo with its ALFA instrument 
llPeek et al.ll200gl , 120091 ; Ivan Loon et~aT]|2009l : iPutman et ail 
2009) are revealing H I structures in our Galactic neighbour- 
hood with unprecedented sensitivity. 

Observational evidence of the interface between atomic 
and molecular gas relies on a good H2 tracer. The J = 
1 — rotational transition of CO is the most common 
surrogate of molecular hydrogen, yet there are conditions 
under whi ch its efficacy as a reli a ble H2 tracer i s un- 
dermined (iDouglas fc Tavlorl [20071: iKlaassen et all 120051 ; 



tional Galactic Plane Survey (IGPS: IStil et all (|2006h : 



iMcClure-Griffiths et all 420051 ): iTavlor et all (|2003l )). Be- 
yond the Galactic Plane, all-sky mapping effo rts of 
the Southern (GASS; IMcClure-Griffiths et al.1 (|2009l )) and 
Northern skies (EBHIS: iKerp et all (| 2009l )) have improved 
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iDobbs fc Bonnelll 120071 : iPelupessv fc Papadopoulod l2009h . 
For example, in photon-dominated regions, the CO molecule 
is largely dissociated, while H2 self-shielding allows for a sig- 
nificant amount of molecular gas to remain. In very cold, 
dense molecular cores, CO will freeze onto dust grains along 
with other larger molecules, and so the H2 content in such 
cores would be underestimated by traditional CO observa- 
tions. In addition to the above scenarios, which assume a 
molecular cloud has already formed, a newly forming molec- 
ular cloud must have an epoch wherein a significant quan- 
tity of H2 has formed before an appreciable (detec table) 
amount of CO has been reached (jBergin et all 12004 ). and 
so the evolution of protostellar objects may be already 
underway in molecular clouds before the onset of criti- 
cal CO densities. OH is another molecule which has been 
used to probe the H-H2 transition region, with some sue- 



2 K. A. Douglas et al. 

cess (|Magnani fc Siskindl [l990l ; IWannier et al.lll993l ). Neu- 
tral and ionised carbon atoms also have p otential to be go od 
tracers of diffuse, pre-CO molecular gas ()Oka et al.l l2005). 

In order for appreciable H2 to form from a largely 
atomic hydrogen reservoir, the gas must be cooled to tem- 
peratures on the order of 20 K, among other favourable 
conditions (sufficient number densities, H2 formati on sites). 
Mode ls of gas flowing into the Perseus spiral arm (|Robertsl 
Il972h predict compression of the H I gas at the interface 
of a "shock ridge." Though his model was isothermal, such 
conditions are favourable for strong cooling of HI, leading 
to self-absorption of 21-cm emission, and possibly enhanced 
H2 production. Hence HI self-absorption (HIS A) observa- 
tions may help to trace the earliest stages of molecular cloud 
formation, in regions where CO is not a workable tracer for 
molecular hydrogen. Widespread HISA h as been observed 
in the Outer Galaxy (jGibson et al.l 120051 ) . and its associa- 
tion with CO is anything but straightforward. HISA and CO 
emission are often found to be coincident spatially and at 
the same radial velocities, but the correspondence is com- 
plicated. There are regions where only one is present, which 
may indicate differences in the evolutionary state of these 
clouds. Studies of HISA and OH emission h ave found the 
OH emission to be more localised to the HISA (|Kavarsll2006l ; 
iLi fc Goldsmithl2003l ). compared to CO. As an extra compli- 
cation, the detection of HISA requires background H I emis- 
sion at the same velocity as the absorbing gas, which may 
not be available along many lines of sight. 

In addition to analytical models, c omputational efforts 
have investigated this scenario in detail. iDobbs et ail (|2008l ) 
performed galaxy-scale SPH simulations that trace the dy- 
namics of interstellar clouds as they encounter shocks in 
spiral arms. To enhance resolution in the spiral arms, their 
model traces only those SPH particles in a galactocentric 
annulus between R g — 5 kpc out to R g — 10 kpc. In addi- 
tion to strong cooling of the atomic gas, molecular hydrogen 
forms in the vicinity of the shocks. Smaller clouds aggregrate 
into larger ones within the spiral shock, where the molecular 
fraction is ^ 0.5. As the clumps leave the spiral arm they are 
sheared into smaller spurs and lose most of their molecular 
gas. 

The link between these models and observations can be 
strengthened by deriving predicted observables directly from 
the simulations. We do not observe densities, temperatures 
and pressures of clouds directly — rather, we rely on spectral 
line and continuum processes to deliver photons containing 
the valuable information, from which we infer the state of 
the object under study. As radiative transfer is a guiding 
principle for this propagation of information, we endeavour 
to simulate the creation of observable photons via radiative 
transfer, within the modelled physi cal representation of the 
ISM in our Galaxy. iGibsonl £2006) demonstrated the pro- 
duction of HISA by spiral shocks with a two-dimensional, 
purely atomic (cold and warm H I) radiative transfer model. 

In this paper we present results from generating syn- 
thetic H I spectral line datacubes from a multiphase numer- 
ic al galaxy calculation . The aforementioned SPH calculation 
bv lDobbs et al.l l|2008l) is used as a basis to create a synthetic 
Galactic Plane survey, fashioned after the H I data products 
created for the IGPS. Simulating such surveys requires the 
"observer" to be placed at an appropriate position inside the 
simulated Galaxy. Such challenges are accomplished using 



TORUS, an Adaptive Mesh Refinement (AMR) code which is 
capable of performing r adiative trans fer calculations to pro- 
duce synthetic spectra l|Harriesll2000l ). Moreover, the added 
flexibility of operating on a numerical model enables us to 
use the same code to observe the Galaxy from multiple po- 
sitions, as well as to emulate extragalactic HI observations 
l|Acreman et al.ll2010l ). 

This paper is organised in the following way. Section 2 
reviews the method used to create our synthetic Galactic 
Plane Survey, and the data cubes produced are described in 
section 3. In section 4 we discuss the use of HISA profiles 
as a diagnostic of early molecular cloud formation. In sec- 
tion 5 we compare briefly our synthetic data to observational 
HI profiles of the Milky Way, setting the basis for a discus- 
sion of how future refinements to our investigation, such as 
the incorporation of feedback mechanisms, will bring closer 
agreement between the synthetic and observed H I distribu- 
tions. 



2 NUMERICAL WORK: SPH SIMULATIONS 
AND TORUS 

The method for cr eating synthetic HI cu bes has been out- 
lined previously in lAcreman et al.l (|2010l ). hereafter called 
Paper I. We briefly summarise the main points here. 

An SPH spiral galaxy calculation is fed into TORUS, 
wherein the particle representation of the gas is transformed 
into a gridded representation, aided by AMR to provide the 
best resolution in regions of highest density. Next, radia- 
tive transfer calculations and ray tracing algorithms produce 
synthetic spectra for lines of sight chosen from an arbitrary 
observing position. The intensities determined by the radia- 
tive transfer and ray tracing are converted to a brightness 
temperature using the Rayleigh- Jeans approximation. In 
Paper I we demonstrated this method by simulating 21-cm 
HI observations of nearby spiral galaxies, assigning the SPH 
galaxies orientations and systemic velocities matching M31 
and M33. 

Moreover, the observer can be placed within the galaxy, 
allowing for the simulation of Galactic HI observations to 
be performed. In such cases, some adaptations of torus 
algorithms must be made. The AMR grid is generated at a 
much finer resolution than that used in Paper I, in order to 
provide better resolution for structure close to the observer. 
The maximum HI mass per cell for the AMR grids used in 
this paper is 2.5 x 10 34 g (12.6 M Q ), which is a factor of 100 
smaller than that used in paper I. In order to cope with the 
increased number of cells in the AMR grid, the domain of 
the grid is restricted to exclude regions of the galaxy which 
are not relevant to the cube being calculated. A different 
AMR grid is used for generating each data cube but they 
typically comprise 6 x 10 unique voxels. 

Also, while the same ray-tracing method is used for 
both the external and internal views, its execution must be 
handled differently if we are to consider multiple lines of 
sight. For the external case the rays are essentially parallel, 
whereas in the internal case there is an obvious angular de- 
pendence to the observations. Additionally, by reversing the 
direction of ray trace, relative to that employed in Paper I 
(i.e. back to front rather than front to back), we are able 
to calculate the contribution of each cell of the AMR grid 
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Figure 1. A synthetic HI spectrum measured from within the 
model galaxy, toward a position near Galactic coordinates I = 
135.5°, b = -0.3°. 



to either the emission or absorption in a given data cube 
pixel. This feature is used to decompose spectral features 
into separate emission and absorption components, as part 
of our investigation of HISA in spiral arms. These differences 
are critical when we extend the creation of single spectra to 
making fully three-dimensional datacubes of simulated HI 
emission. The cubes are discussed in the next section. 

Figure [T] shows an H I spectrum measured from a po- 
sition within the SPH galaxy, corresponding roughly to the 
Earth's position in our Galaxy, toward a direction corre- 
sponding to Galactic longitude I » 135.5° and latitude 
b ~ —0.3°. We see "local" spiral arm gas in a radial veloc- 
ity range —20 < v r < km s^ 1 , and an outer arm beyond 
v r ~ —38 km s _1 . The outer arm shows some evidence of H I 
absorption, with a narrow dip in XL appearing near v r = —41 
km s _1 . This HISA signature is common in synthetic spec- 
tra all across our region of interest. Since HISA is a main 
focus of our investigation into molecular cloud formation, 
the TORUS code was modified so that in addition to total 
brightness temperature cubes, we create separate datacubes 
where only positive or negative intensity contributions (dl) 
for a given grid cell are mapped. 

For this paper, the simulated galaxy has reached a time 
of 250 Myr. At this age the distribution of gas in different 
phases, and the fraction of molecular gas, have reached a 
roughly steady state. While tracing the evolution of specific 
ISM structures using many timesteps from an SPH simula- 
tion is the subject of future investigations, the present work 
concentrates on deriving HI properties from a single epoch, 
as is the case with most Galactic radio observations. As men- 
tioned previously, the models are restricted to Galactocen- 
tric radii of 5 to 10 kpc. For this reason it is reasonable to 
aim to reproduce observations of the Outer Galaxy, toward 
its Anticenter (I — 180°). In this paper, we survey the en- 
tire second quadrant of our SPH galaxy, covering longitudes 
between I = 90° and i = 180°. 
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Figure 2. A midplane density slice of the simulated Galaxy grid- 
ded by TORUS, and showing the position of our "observer" for the 
synthetic Galactic Plane Survey. 



Table 1. Second Quadrant Mosaic Regions 



Region Name 


^■central ( ) 


ABC 


99.0 


DEF 


117.0 


GHI 


135.0 


JKL 


153.0 


MNO 


171.0 



3 SYNTHETIC DATACUBES OF H I 
BRIGHTNESS 

For the Galactic point of view, we construct Tb(£, b, v r ) cubes 
with angular coordinates corresponding to "Galactic" coor- 
dinates from the position of the observer. This position is 
chosen to be in the midplane, 2.2 x 10 22 cm (7.13 kpc) from 
the galaxy's center, as illustrated in figure[2] This choice cor- 
responds to a position in an interarm "spur" of H I, similar 
to the Orion Spur of our Milky Way. 

3.1 Channel Maps 

A total of fifteen cubes of angular dimension 6.08° x 6.08° are 
produced to cover the entire Second Quadrant, with a small 
amount of overlap between adjacent fields for mosaicking 
purposes. With an image scale of 0.5' per pixel, each channel 
map produced by TORUS measures 732 2 pixels. Our spectral 
range covers radial velocities from +30 km s _1 to —120 km 
s" 1 sampled at 0.5 km s -1 , resulting in 300 channels for 
each cube produced by TORUS. 

For ease of visualisation, the fifteen cubes are combined 
into five mosaics, summarised in table [T] Each mosaic re- 
gion therefore covers just over 18° of Galactic longitude, 
centered on the value of Icmtrai given in the table. For il- 
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Figure 3. Synthetic Perseus Arm (SPA) gas near v r = —45 km/s 
for the GHI region (126° s$«s$ 144°). 



lustrative purposes we show a channel map from the central 
mosaic, region GHI, in Figure [3] This channel map shows 
HI gas near v r = —45 km s _1 , as our aim is to study HISA 
produced in spiral arm gas with radial velocities correspond- 
ing to Perseus Arm gas (we will refer to the simulation's 
"Perseus Arm" as the "SPA"). 

A common feature seen in our SPA channel maps is 
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Figure 4. Column density map (top) and optical depth channel 
map (bottom) for the SPA cube centered on I = 135° ,6 = 0°. 



that the brightest HI emission appears not in the central 
midplane (b — 0°) of the spiral arm, where the densities 
are greatest, but in ridges at approximately ±1° above and 
below the midplane. From the known distribution of gas col- 
umn density in the original SPH simulations, the expected 
HI distribution assuming optical thinness would not have 
predicted such a pattern for maximal brightness temper- 
atures. Rather, this is an absorption effect arising from a 
highly dense and cold central (|6| < 0.5°) midplane. As fig- 
ure^ (top panel) demonstrates for the cube toward £ — 135°, 
the synthetic HI clouds occupy the plane quite randomly, 
mostly within |6| < 1°. With an increase in column density, 
an increased optical depth r of H I follows naturally. As a re- 
sult, the midplane of the SPA is largely optically thick. With 
TORUS we can construct cubes of r(£, b, v r ) in addition to Tt, 
and we show an example channel map of r in the SPA for 
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the same region in the lower panel of figure [4] Individual HI 
clouds exhibit optical depth values r 15, and thus the HI 
gas nearer the midplane absorbs much more strongly, and 
the brightness temperatures are diminished by the combi- 
nation of high opacity and colder temperatures. The optical 
depths above |b| = 0.5° are considerably lower, and so we 
observe emission rather than absorption. 

Another observed (and related) consequence of the over- 
dense, high-r SPA plane is that the amount of absorption 
seen in the SPA greatly exceeds the HISA distribution found 
in the Perseus Arm of the Milky Way. While this level of self- 
absorpti on is greater than what is expected from Galactic 
studies (|Gibson et al.ll2005l ; [Goldsmith et al.ll2007l ). the ex- 
aggerated contrast between emission and absorption in the 
SPH galaxy allows us to investigate the HISA phenomenon 
in a more "controlled" interstellar environment. As men- 
tioned above, we make separate cubes which track pixels 
where brightness temperature increases or decreases along 
the line of sight, as a means of tracing regions where self- 
absorption occurs. 

Figure O shows positive-only and negative-only spectra 
for the line of sight total- Tf, spectrum shown in figure [T] and 
channel maps for region GHI at the same radial velocity as 
shown in figure [3] The negative-only spectra and channel 
maps are very useful in isolating HISA features in our syn- 
thetic survey; the majority of negative- Tt pixels lie in the 
plane at radial velocities corresponding to the SPA. We will 
return to our analysis of HISA features in the next section. 
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3.2 Position- Velocity Views 

In figure[6]we show the Galactic view toward t — 135° again 
(region GHI), with an I — v r cut along b — 0°. The main 
feature to be elucidated is the separation of local (v r w 0) 
HI from gas at synthetic Perseus Arm (SPA) velocities, near 
—35 to —55 km/s. The local gas has little discernable struc- 
ture in our synthetic cube, as its proximity to the observer 
position means the gas has a very large angular extent. As 
is the case with our Galaxy, the mean radial velocity of the 
SPA becomes less negative toward higher longitudes, and 
blends with the local gas toward the Anticenter, a projection 
effect based on Galactic rotation. Toward lower longitudes, 
the SPA spans a larger range of radial velocities. 

Beyond the SPA, there is no outer spiral arm across 
the Second Quadrant, until the appearance of terminating 
spiral arms from spiral structure originating toward the In- 
ner Galaxy. Across the Second Quadrant there are discrete 
clouds, containing high fractions of molecular gas in many 
instances, which exit the outer edge of the SPA, but these 
clouds mostly share radial velocities with H I gas within the 
SPA. 

As we move away from the midplane in the longitude- 
velocity projection, the SPA shows bright HI emission out 
to |6| « 1.3°, after which its intensity drops sharply as one 
moves further from the midplane. This effect is shown to- 
ward t = 135° in figure [3 With a distance to the SPA of 2 
kpc (at its closest), this vertical distribution suggests that 
most of the HI is confined to within ~ 100 pc of the mid- 
plane. 
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Figure 7. Latitude- velocity slice of synthetic H I toward a Galac- 
tic longitude of i = 135° . The SPA H I emission is mostly confined 
to within 1 degree of the midplane. 



4 HISA AND SPIRAL SHOCKS 

In our Galaxy, HISA features tend to outline individ- 
ual clouds or complexes of cold H I ([Gibson et al. 1 120051 ; 
iKavars et all 120051 ; llvIcClure-Griffiths et alj 120061 ). By con- 
trast, our synthetic data show large bands of HISA that 
subtend much larger angular extents. While unrealistic, this 
large abundance of HISA gives us many sightlines to test its 
origin. 
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Figure 5. Synthetic HI emission-only (left) and absorption-only (right) spectra and channel maps. The spectra are for the same line of 
sight as in figure [l] but with the positive and negative 7], contributions separated. The channel maps show the same, for all of region 
GHI (see figure [3}. 
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Figure 6. Longitude- velocity diagram of simulated HI emission near i = 135°, b = 0° for the Galactic synthetic observation. Local and 
Perseus Arm gas are well separated. 




Distance from observer (kpc) 

Figure 8. Radial velocity of the gas plotted against distance, for 
I K, 130.5°, b w -0.38° in the simulated galaxy. 



The scenario described by iRobertsI j 19721) for shock- 
induced compression of H I (as evidenced via HIS A) can be 
tested directly with our simulations. Figure [8] shows a plot 
of the radial velocity for gas along a line of sight toward a 
position in region GHI, with I ~ 130.5° and 6 w —0.38°. 
Between the observer and a distance of about 2 kpc, the 
radial velocity of the gas decreases steadily as a result of 
the gas moving with the rotation of the galaxy. Then, at 
a distance of 2 kpc, the radial velocity undergoes a sharp 
decrease from about —40 km s _1 down to a local minimum 
of approximately —50 km s _1 . The radial velocity remains 
in this range as we continue to move away from the observer, 
until about 2.2 kpc, where it slowly rises to —40 km s" 1 and 
levels off for about half a kiloparsec, before again decreasing 
very gradually beyond 2.7 kpc. 

The dip near 2 kpc corresponds to the near edge of the 
SPA, and we interpret this sharp decline in radial velocity 
as the result of the H I gas experiencing a shock in the spi- 
ral arm, where it loses much of its line of sight velocity as 
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Figure 9. The HISA profile (ie. negative- Tf,-only spectrum) to- 
ward the position described in figurc[8] The radial velocities where 
HISA occurs are between —40 and —50 km s _1 , exactly where it 
would be expected from figure [8] 



it begins to travel along the spiral arm edge. The gas cools 
during this process, and so the hydrogen at 2 kpc will absorb 
the H I emission from the gas behind it with the same radial 
velocity, resulting in strong HISA present in the spectrum 
between —40 and —50 km s . Indeed we show in figure [5] 
that this is the precise radial velocity range where HISA is 
prevalent at thi s same position in the datacube. This con- 
firmation of the lRoberts! (|l972T ) scenario for HI shocks, in a 
region where molecular clouds are also formed in the model, 
shows the value of HISA as a tracer of the atomic-molecular 
interface of the ISM. 

In figure [TO] we zoom in to distances corresponding to 
the SPA, for the same line of sight as figure [8] While this 
sightline shows appreciable HISA, there is nothing extraor- 
dinary about it compared to other regions across the SPA. 
It is a fair representation of the type of H I profiles observed 
across the midplane of our synthetic survey. As shown in the 
top left panel, the density of neutral hydrogen rises by two 
orders of magnitude between 2.0 and 2.1 kpc from the ob- 
server, as the spiral arm boundary is crossed. While a single- 
fluid (7 = I ) adiabatic shock would result in at most a factor 
of four density enhancement, such conditions would not con- 
strain a strongly cooling shock, which describes the present 
scenario. Indeed, a corresponding drop in gas temperature is 
shown in the top right panel, as the gas cools strongly from 
values of ~ 3000 K down to less than 100 K. In the lower 
left figure, we plot the H2 fraction, by mass, of the total 
gas. The first peak, between 2.0 and 2.1 kpc, is interpreted 
as newly formed H2 that has been formed in the vicinity of 
the shock at that site in the SPA. The other two high-H2 
areas beyond the shock, near 2.25 and 2.6 kpc, would be 
(giant) molecular clouds that were formed upstream by the 
same shock mechanism, at higher longitude, and which are 
moving through the SPA. In the absence of feedback in our 
simulation, these clouds may be larger than expected, and 
longer-lived. The outer edge of the SPA is less easy to de- 
fine than the inner edge, given that there is no sharp outer 
boundary, but across much of the SPA, the density begins 
to decrease toward interarm values approximately 0.7 to 0.8 
kpc behind the inner boundary. 



The effect of these density and temperature variations 
on the observed H I profile can be informed by the lower right 
panel of figure [TCT1 which shows the intensity contribution dl 
per unit of column density as a function of distance. HISA 
(negative dl) appears to be produced behind the shock, once 
the temperature has dropped to below 100 K, and where 
the HI density is above 10 -23 g cm -3 . At the front of the 
shock there is a narrow region where emission is produced 
rather than absorption, because the gas is still quite warm, 
yet the HI density is already increasing rapidly. Such con- 
ditions are favourable for the onset of H2 production, which 
depen ds strongly on den sity as well as weakly on temper- 
ature jBergin et all 12004") ■ Still, this shows that within the 
region where the gas velocity is strongly influenced by the 
spiral arm shock, both emission and absorption from HI 
atoms are possible. Yet we can constrain the occurrence of 
HISA to cold, dense regions which also contain molecular 
gas. The temperature appears to remain at values around 
100 K out to a distance of ~ 2.65 kpc, after which it begins 
to rise and the molecular fraction drops sharply. Thus the 
cold gas is mostly confined to the spiral arms, a nd th e gas is 
warmed in the interarm regions. iDobbs et alj l|2008l ) found 
that roughly 62 per cent of the H I was in the cold phase, 
where T ^ 150 K was their criterion for cold gas. 

The shock region itself is much smaller than the ~ 100 
pc region behind it where the physical state of the ISM un- 
dergoes these interesting changes. In the original SPH sim- 
ulations, the smoothing length decreases greatly in high- 
density regions, dropping to values of about 5 to 6 pc to 
resolve region s where the density approaches 10~ 22 g cm -3 
|Dobbsl 120081 ). which is approximately the peak density 
achieved for the shocked HI shown in figure [TOl SPH simula- 
tions typically require several smoothing lengths to capture 
a shock, and the shock width in any numerical mod el will 
be co nsiderably greater than any realistic shock width (|Pricel 
l2008h . 



5 DISCUSSION 

5.1 Comparison with Observational Data 

The synthetic datacubes we have produced correspond to 
the Second Galactic Quadrant, for which the best compar- 
ison set is the Canadian Galactic Pl ane Survey, part o f the 
International Galactic Plane Survey (|Tavlor et al.ll2003T ). We 
compare below the general characteristics and features ob- 
served in both the synthetic and real data. By the very fact 
that our SPH galaxy is not the Milky Way, no strong agree- 
ment is expected a priori. The main purpose of the synthetic 
models is to investigate whether the geometry, kinematics, 
chemistry and thermodynamics of a spiral shock can pro- 
duce HISA, as has been observed in the Perseus Arm of 
our Galaxy. In comparing the two datasets, we can discern 
what characteristics this relatively simple galaxy model has 
in common with Galactic HI observations. 

In figure QT] we show the Galactic analogue of the syn- 
thetic H I channel map for region GHI that was shown in fig- 
ure(3] The main difference observed across the Second Quad- 
rant is the apparent confinement of H I gas to the plane in 
our synthetic data cubes, compared to the widespread and 
more diffuse distribution of atomic gas seen in the CGPS 
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Figure 10. A close-up of the shock region seen in figure[8] The shock front position is approximated by the vertical dashed line. As HI 
density increases (top left) and temperature decreases (top right) , H2 is formed in the shock region (bottom left), and background HI 
emission is absorbed (negative dl values in bottom right panel). 



data. This is not surprising, given the lack of feedback mech- 
anisms in the simulations that would serve to stir up the 
interstellar medium more vigourously — we will discuss this 
point further below. For now, we consider the consequences 
of this scale-height mismatch. For one, with a greater column 
density of HI toward the overdense midplane, the bright- 
ness temperatures for the synthetic line profiles are larger 
than what is observed in the Milky Way. Within the CGPS 
data brightness temperatures above 150 K are not common, 
yet in the synthetic cubes we encounter many regions with 
T b > 200 K. 

Unlike the scale height differences observed in the chan- 
nel maps, the spiral structure in the SPH galaxy and the 
Milky Way show good agreement, in terms of the Perseus 
Arm, the main feature in the Second Quadrant of the Outer 
Galaxy. Figure [12] shows the Milky Way's mid-Second Quad- 
rant midplane longitude-velocity distribution, which por- 
trays a similar separation, on the order of Av r ~ 40 km 
s _1 , between local and Perseus Arm gas to that observed in 
figure [6] The presence of gas beyond the Perseus Arm is a 
main difference exhibited by the real data, as well as distinct 
HI features in the interarm gas. Additionally, the overall ve- 
locity width of the Galaxy's spiral arms is in general greater 
than that of the simulated galaxy. In the SPH galaxy, parti- 
cles were given circular velocities with a velocity dispersion 



assigned to be a = 6 km s 1 , while values of a ~ 8 km s 1 
or greater typica lly describe the HI of our Outer Galaxy 
jSaha et alj"l2009l V Many factors contribute to such a dis- 
persi on value, includ i ng fee dback mechanisms and a warped 
disk; iDobbs fc Price! (2008) also found that magnetic fields 
broaden spiral arms in SPH galaxy simulations. 



5.2 The Role of Feedback 

The interstellar medium of our Galaxy is constantly dis- 
rupted, by both external and internal agents. The inflow 
of gas from satellite galaxies and high-velocity clouds pro- 
vides fuel for future generations of molecular clouds and 
stellar systems. Numerous feedback processes from stellar 
evolution also add energy and momentum to the ISM from 
within: outflows, supernova shocks and winds from individ- 
ual stars or stellar associations. With these agents absent 
from the simulated galaxy, our synthetic observations cor- 
respond to a rather quiescent version of a Milky Way-type 
galaxy. The consequences are abundantly clear in our syn- 
thetic datacubes. HI clouds collect toward the midplane and 
cool, with nothing to destroy them, resulting in larger than 
expected column and number densities. Thus the lower scale 
height seen in our simulated galaxy is well explained. If 
our synthetic data had shown agreement with the HI scale 
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Figure 12. Longitude- velocity diagram of Milky Way HI emission near t = 135°, 6 = 0° as observed with CGPS data. 



height suggested by the CGPS data, then the importance of 
feedback in contributing to the overall structure of the ISM 
would be severely diminished. 

Much of the HI distribution in our Galaxy can be de- 
scribed with words such as filaments, shells, and bubbles. 
Interstellar bubbles are created by outward forces originat- 
ing from a central nidus. Shells and filaments form as a result 
of bubbles encountering each other, or other ISM structures, 
and their evolution is also aided by radiative and magnetic 
field effects. Without the injection of supernova-like events 
or winds from OB associations, the simulated galaxy cannot 
be expected to recreate the filamentary structure observed 
with H I observations of the Milky Way. Hence feedback is a 
possible solution to both the overconfinement of H I in the 
midplane, as well as the non-filamentary nature of the HI 
channel maps. Feedback should also increase the velocity 
dispersion of the simulated galaxy's spiral arms, since the 
mixing of the HI will certainly result in gas with a wider 
range of radial velocities. 

5.3 Future Work 

The simulation used in this paper provides a basis with 
which to compare future simulations which include more 
physical processes. As highlighted above, feedback appears 
to be instrumental in reproducing the correct structure 
(holes, filaments) and vertical distribution of the ISM. More- 
over, magnetic fields are an important sour ce of press ure, 
supporting the gas in the vertical direction !|Coxll2005D . In 
addition, magnetic fields are l ikely to reduce thermal in- 
stability in the gas (|Fieldlll965l ; iHennebelle fc Peraultlbood; 
lHeitsch"etT al, 2009) , which could reduce the amount of HISA 
seen in our simulated observations. On the other hand, self 
gravity will promote the formation of dense HI or molecu- 
lar structures in the midplane, particularly at higher surface 
densities. Our synthetic observations of the HI distribution 
and the resulting HISA properties in these future simula- 
tions will demonstrate the cumulative effect of the added 
physics. 

As the molecular fraction of the SPH galaxy evolves, 
formation of H2 and CO tak e place using prescribed rate 
equations (|Dobbs et alj |2008). The simplicity of the two- 
level hyperfine 21-cm transition makes the production of 
the synthetic HI spectra rather straightforward, compared 
to the detailed balance calculations required for other ISM 
tracers. Still, in principle the TORUS code can be adapted to 
produce synthetic o bservations of any sp ecies present in the 
original simulation. iRundle et al.l (|2010l ) have achieved syn- 



thetic molecular line observations for several species, includ- 
ing CO, for far-field view simulated observations of molecu- 
lar clouds with TORUS. We will integrate their work to pro- 
duce a synthetic CO galactic plane survey to complement 
our H I data. Besides spectral lines, we will also aim to pro- 
duce realistic maps of continuum emission from dust. 

The use of numerical simulations also offers the unique 
opportunity to trace individual feat ures in the time do main. 
The simulation we used here from iDobbs et all |2008) and 
others like it show the evolution of atomic hydrogen clouds 
as they encounter shocks in spiral arms, form molecular 
clouds, and eventually disperse as they exit the spiral arm 
on the other side. Radio observations can only piece together 
this scenario by observing multiple clouds at different stages 
of this evolution. In a future study we will recreate our syn- 
thetic Galactic Plane survey for several timesteps of a nu- 
merical solution, and track the evolution of the HI profiles 
of gas clouds which pass through spiral shocks. 

In all of the above studies, we can also extend our syn- 
thetic Galactic Plane survey to other parts of the simulated 
Galaxy. From our chosen observer position, only lines of 
sight within \l\ < 46° are affected by the absence in the 
model domain of the inner 5 kpc of the galaxy. Thus we can 
extend our investigation to include H I observations from 
other components of the IGPS for comparison. 



6 SUMMARY 

We have demonstrated our ability to conduct a synthetic 
Galactic Plane HI survey. With the aid of the versatile 
TORUS code, we converted an SPH simulation of a spiral 
galaxy section into a set of HI datacubes predicting the 
brightness temperature distribution from an observer's posi- 
tion within that galaxy. With these cubes we can trace global 
features such as HISA and spiral arm structure, and in com- 
parison with Galactic Plane Survey cubes of our Galaxy, we 
find common features, as well as clear and explicable dif- 
ferences in the HI distribution. The main differences seen, 
which include a lack of filamentary structure, a lower HI 
disc scale height, and spiral arms with a lower velocity dis- 
persion, may be well explained by the absence of feedback 
processes that would conspire to mix the SPH galaxy's ISM 
more vigourously. We have shown that HISA is produced 
naturally by the influence of a spiral arm in the outer galaxy, 
where little or no HISA would otherwise be expected, be- 
cause the unperturbed line-of-sight velocity field is mono- 
tonic with distance. With future SPH simulations that in- 
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Figure 11. CGPS Channel map of Perseus Arm gas centered on 
t = 135°, b = 0°. 



elude feedback mechanisms, we will discern whether stellar 
processing of interstellar gas leads to a closer match between 
the model and actual galaxy. We will also track the evolu- 
tionary trajectory of atomic clouds seen as HISA, to inves- 
tigate further the early stage of molecular cloud formation. 
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